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ABSTRACT 

O ■ By definition, a Type Ic supernova (SN Ic) does not have conspicuous lines 

^ ' of hydrogen or helium in its optical spectrum. SNe Ic usually are modelled in 

J>^' terms of the gravitational collapse of bare carbon-oxygen cores. We consider the 

possibility that the spectra of ordinary (SN 19941-like) SNe Ic have been misinter- 
preted, and that SNe Ic eject hydrogen. An absorption feature usually attributed 
to a blend of Si II A6355 and C II A6580 may be produced by Ha. If SN 19941-like 
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CN ■ SNe Ic eject hydrogen, the possibility that hypernova (SN 1998bw-like) SNe Ic, 

(-^ ■ some of which are associated with gamma-ray bursts, also eject hydrogen should 

^ ■ be considered. The implications of hydrogen for SN Ic progenitors and explosion 

Tj- ■ models are briefiy discussed. 

O" 

O ■ Subject headings: supernovae: general — supernovae: individual (SN 19941, SN 

1999ex) 

O 

U 

^ ■ 1. INTRODUCTION 



X 



A Type II supernova (SN II) has conspicuous hydrogen fines in its optical spectrum. A 
Type lb supernova lacks conspicuous fiydrogen lines but does fiave conspicuous He I lines. 
A Type Ic supernova fias conspicuous lines of neitfier fiydrogen nor He I. For a review of 
supernova spectral classification see Filippenko (1997). 

Supernovae of tfiese types are tfiougfit to result from core-collapse in massive stars. A 
SN lb is fiydrogen-deficient in its outer layers, and displays He I lines owing to nontfiermal 
excitation by tfie decay products of radioactive ^^Ni and ^^Co (Lucy 1991). In recent years 
it fias becomes clear tfiat at least some SNe lb are not completely fiydrogen-free; in fact, 
it is likely tfiat most or even all SNe lb eject a small amount of fiydrogen at fiigfi velocities 
(Deng et al. 2000; Brancfi et al. 2002; Elmfiamdi et al. 2006). If some of tfie events 
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classified as Type lb had been observed earlier, hydrogen hnes might have been conspicuous, 
in which case they would be classified as Type lib — the designation used for events such as 
SN 1993J that looked like a Type II at early times but later looked like a Type lb (Filippenko, 
Matheson, & Ho 1993). 

A SN Ic either is helium-deficient in its outer layers or fails to nonthermally excite its 
helium (Woosley & Eastman 1997). It has become common practice to model SNe Ic in 
terms of core collapse in bare (or nearly bare) carbon-oxygen cores (Iwamoto et al. 1994; 
Foley et al. 2003; Mazzali et al. 2004). 

Supernova spectral features are P-Cygni features characterized by line-centered emis- 
sion components and blueshifted absorption components, with the absorptions frequently 
being more identifiable. Spectral features usually are blended owing to the huge Doppler 
broadening. Nevertheless, some identifications of spectral lines are definite. For example, 
blends of Fe II lines and features produced by Ca II H&K and the Ca II infrared triplet 
appear in all types of supernovae as long as the temperature is sufficiently low. However, 
there also are some serious identification ambiguities. The most well known involves the 
Na I D-line doublet at mean wavelength A5892 and the strongest optical line of He I, A5876, 
which are separated by only about 800 km s^^. If the corresponding observed feature is 
strong and other He I lines are not present, the feature is produced at least mainly by Na I, 
while if other He I lines are clearly present then the feature is at least partly due to He I. 
But if the observed feature is not strong, it can be difficult to choose between Na I and He I. 

When SNe Ic are interpreted in terms of bare carbon-oxygen cores, an absorption feature 
usually near 6200 A, which we will refer to as the 6200 A absorption, is attributed to the 
strongest optical line of Si II, A6355 (the transition definitely responsible for the similarly 
located deep absorption in Type la supernovae), perhaps blended with the strongest optical 
line of C II, A6580, forming in higher-velocity ejecta than Si II. Local-thermodynamic- 
equilibrium (LTE) calculations of Sobolev line optical depths for a composition dominated 
by carbon and oxygen (with hydrogen and helium burned to carbon and oxygen, and solar 
mass fractions of heavier elements) show that within certain intervals of temperature Si II 
A6355 and C II A6580 are expected to have significant optical depths (Hatano et al. 1999). 
Nevertheless, the identification of the 6200 A absorption is plagued by ambiguities. The 
strongest line of Ne I, A6402, is about 2200 km s~^ to the red of Si II A6355, and Ha 
A6563 is only about 800 km s^^ to the blue of C II A6580. Because these four ions, each of 
which could appear in supernova spectra (although Ne I probably would require nonthermal 
excitation), have their strongest optical lines to the red, but not too far to the red, of the 
6200 A absorption, the ambiguity is difficult to resolve when the observed feature is weak 
and other lines of these ions do not produce identifiable features. 
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In this paper we are primarily concerned with the possible presence of hydrogen in 
SNe Ic. This issue (and that of helium in SNe Ic) has been addressed in the literature 
before; for summaries see Filippenko (1997); Matheson et al. (2001); and Branch (2002). 
The presence of hydrogen has been suggested, for example, by Filippenko (1988, 1992); 
Filippenko, Porter, & Sargent (1990); Jeffery et al. (1991); and Branch (2002). The currently 
prevailing view, however, is that hydrogen is absent and the 6200 A absorption is produced 
by Si II and/or C II (Wheeler et al. 1994; Millard et al. 1999). In §2, based on a comparison 
of spectra of the Type lb or Type Ib/c SN 1999ex, which we believe to contain hydrogen, and 
the Type Ic SN 19941, we raise the question of whether hydrogen also is present in SN 19941. 
In §3, comparisons of spectra of SN 1999ex and SN 19941 with synthetic spectra generated 
with the parameterized resonant-scattering code Synow are presented and discussed. In §4, 
the implications of the presence of hydrogen in ordinary (SN 19941-like), and possibly even 
in hypernova (SN 1998bw-like; Foley et al. 2003) SNe Ic, are briefly considered. 



2. COMPARISON OF SN 1999ex AND SN 19941 

SN 1999ex was initially classified as Type Ic based on a resemblance of its early spectrum 
to SN 19941 (Hamuy & Phillips 1999), but because it later developed He I lines Hamuy et al. 
(2002) revised the classification to intermediate Type Ib/c. Although the He I lines were 
not as deep as in most events that have been classified as Type lb, their presence was 
definite, so Branch (2002) referred to SN 1999ex as a "shallow-helium" Type lb. Hamuy 
et al. (2002) labelled the 6200 A absorption as Si II, but Branch (2002) argued that the 
correct identification is Ha, consistent with the presence of Ha in most if not all other 
SNe lb. Recently Elmhamdi et al. (2006) have reinforced the conclusion that SNe lb, 
including SN 1999ex and another transition Type Ib/c event, SN 1996aq, eject hydrogen. 
(The overluminous Type lb SN 1991D (Benetti et al. 2002) could be an exception.) 

The spectra of SN 1999ex and the Type Ic SN 19941 shown in Figures 1 to 3 are 
from Hamuy et al. (2002) and Filippenko et al. (1995), respectively. The horizontal scale 
on all figures in this paper is logarithmic wavelength, which allows the widths of Doppler- 
broadened features to be compared on an equal basis across the whole spectrum. The spectra 
are "tilted" by multiplying by A", with a chosen to make the flux peaks near 4600 A and 
6300 A about equally high. The tilting makes it easier to compare spectral features. When 
comparing SN 1999ex and SN 19941 spectra we artificially blueshift those of SN 1999ex in 
order to compensate for the different photospheric velocities and roughly align the absorption 
features. 

In Figure 1 a spectrum of SN 19941 obtained 2 days before the time of maximum light 
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in the B band (day —2) is compared with a day —1 spectrum of SN 1999ex that has been 
blueshifted by 4000 km s~^. Figure 2 is hke Figure 1, but for day +4 spectra of SN 19941 and 
SN 1999ex, with the latter blueshifted by 2000 km s~^. Figure 3 is for a day +10 spectrum of 
SN 19941 and a day +13 spectrum of SN 1999ex, with the latter blueshifted by 2000 km s^^. 
In all three figures the spectra are similar in many respects, although the features are more 
washed out in SN 19941 because of its higher photospheric velocity at each epoch (hence the 
necessity to blueshift SN 1999ex to illustrate the similarities). Otherwise the main differences 
are that He I A6678 and A7065, although weak, are clearly present in SN 1999ex (making it 
a Type lb or at least a Type Ib/c), but they are not clearly present in SN 19941 (making it 
a Type Ic). However, it is not possible on the basis of Figures 1 to 3 to exclude the presence 
of weak He I features in SN 19941. Given the similarities in these figures — and they are 
rather striking — the question becomes: are the Na I and Si H/C H identifications correct 
for SN 19941? If so, then the resemblance of the SN 1999ex and SN 19941 spectra from about 
5500 A to 6300 A is coincidental. 



3. COMPARISONS WITH SYNTHETIC SPECTRA 

To further investigate the possibility of hydrogen (and He I) in SN 19941 we have used the 
parameterized resonant-scattering Synow code (Branch et al. 2002) to generate synthetic 
spectra for comparison with spectra of SN 1999ex and SN 19941. 



3.1. SN 1999ex 

In Figure 4 the day +4 spectrum of SN 1999ex is compared with a synthetic spectrum 
in which the 6200 A absorption is produced by Ha. In the synthetic spectrum the velocity 
at the photosphere is 8000 km s~^. Hydrogen is detached from the photosphere at 13,000 
km s~^ where Ha has an optical depth of 0.5. Hydrogen-line optical depths decrease outward 
exponentially but very slowly, with e-folding velocity Ve = 20,000 km s^^, to an imposed 
maximum velocity of 20,000 km s^^. Given the closeness of the fit to the 6200 A absorption, 
we know that with fine adjustments of the optical-depth profile we could refine the fit to 
be practically perfect. Thus at the Synow level of analysis. Ha is completely adequate and 
we would not prefer a different identification without other evidence. For all ions used in 
Figure 4, reference-line optical depths and Ve values (in units of 1000 km s~^) are given in 
Table 1. 

Elmhamdi et al. (2006) show that Si II A6355 cannot account for the 6200 A absorption 
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of SN 1999ex on its own because even when it is undetached from the photosphere, its 
synthetic absorption feature is too blue. It is unhkely that the 6200 A absorption is due 
mainly to Ne I A6402 because the observed feature is strong enough that when it is fitted 
with Ne I A6402 other Ne I lines appear in the synthetic spectrum and the fit deteriorates. 
This conclusion is based on LTE excitation and must be tested against detailed non-LTE 
calculations that take nonthermal excitation into account. 

Hamuy et al. (2002) also obtained a day +4 spectrum of SN 1999ex that extends 
to nearly 2.5 microns. Synow fitting parameters used to fit the optical He I lines also 
give satisfactory fits to features attributed to He I A10830 and He I A20851 (Branch 2002; 
Elmhamdi et al. 2006), but the infrared spectrum does not provide a useful constraint on 
the presence of hydrogen. 



3.2. SN 19941 

Synow fits to spectra of SN 19941 in which the 6200 A absorption is attributed to 
Si H and/or detached high-velocity C II are shown in Millard et al. (1999) and Elmhamdi 
et al. (2006). For this paper we began trying to fit the day +4 spectrum of SN 19941 by 
varying the input parameters that were used for the synthetic spectrum shown in Figure 4 for 
SN 1999ex, i.e., we assumed the presence of hydrogen and He I lines. Figure 5 for SN 19941 
is like Figure 4 for SN 1999ex. The synthetic spectrum has a velocity at the photosphere of 
11,000 km s~^ (compared with 8000 km s^^ for SN 1999ex). Hydrogen lines are detached 
from the photosphere at 15,000 km s^^ where Ha has an optical depth of 0.4. Hydrogen-line 
optical depths decrease outward exponentially with e-folding velocity Ve = 20,000 km s~^ 
to an imposed maximum velocity of 22,000 km s~^. This hydrogen optical-depth profile is 
quite similar to that used for SN 1999ex in Figure 4, and in Figure 5 the 6200 A absorption is 
fit nicely. The 5700 A absorption also is fit nicely by He I A5876. In the synthetic spectrum 
He I A6678 and A7065 are weak and rather washed out owing to the high photospheric 
velocity, but they do more good than harm and their presence in the observed spectrum 
cannot be excluded. (However, regarding the strong observed absorption near one micron, 
the situation remains as illustrated and discussed in Millard et al. (1999); the synthetic He I 
A10830 absorption is too weak and narrow to account entirely for the observed absorption.) 
Reference-line optical depths and Ve values for Figure 5 are given in Table 1. 

The synthetic spectrum in Figure 6 is like that of Figure 5 except that it includes 
Si II lines instead of hydrogen lines. The figure illustrates that in SN 19941, as well as in 
SN 1999ex, Si II A6355 is too blue to account for the 6200 A absorption on its own (fiux 
differences between observed and Synow spectra are to be expected, but differences in wave- 
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lengths of absorption features are regarded as serious disrepancies) , although it is possible 
to obtain a fit by invoking detached high-velocity C II. Because the 6200 A absorption is 
less deep in SN 19941 than in SN 1999ex, and the SN 19941 features are more smeared out 
by the higher photospheric velocity, it also is difficult to exclude Ne I A6402 in SN 19941. 



4. DISCUSSION 

4.1. Spectroscopy 

We have discussed the possibility of Ha in absorption, in the spectrum of the Type Ic 
SN 19941. As shown in Figure 7, this is an issue for SNe Ic in general; the spectra within the 
wavelength interval emphasized in this paper have strong similarities. Other evidence for 
hydrogen in SNe Ic was emphasized by Filippenko (1988^; 1992): the apparent presence of 
broad Ha emission in near-maximum-light spectra. This feature can be seen in the day —2 
spectrum of SN 19941 in Figure 1 (beneath a narrow Ha emission due to an H II region). In 
order to produce a peaked rather than a flat Ha emission component the hydrogen would 
have to be present down to the photosphere, rather than being confined to a detached high- 
velocity shell such as we have invoked for the absorption component. In this regard it is 
interesting to note that in Figures 1 to 3 the red edges of the putative Ha and He I A5876 
absorptions in SN 19941 are less sharp than in SN 1999ex; a detached shell produces sharp 
red edges. A continuous Ha optical-depth distribution that peaks at 15,000 km s~^ but 
extends down to the photosphere might suffice to produce a peaked Ha emission. This can 
only be tested by non-LTE calculations because the resonance-scattering approximation on 
which Synow is based is known to be a poor approximation for Ha in SNe II. It is not 
clear that a large hydrogen mass would be required because the appearance of hydrogen 
and He I lines depends on the radial optical-depth profiles (Branch et al. 2002) and on the 
distribution of ^^Ni (Woosley & Eastman 1997). 

At present it is difficult to establish the presence or absence of hydrogen in SN 19941- 
like SNe Ic, but the issue should be considered further. A large-scale comparative study 
of observed spectra, like that of Elmhamdi et al. (2006) but including SNe Ic as well as 
SNe lb, probably would yield clues. In addition, detailed non-LTE spectrum calculations 
for a grid of explosion models are needed. Such calculations, even based on parameterized 
explosion models rather than full hydrodynamical calculations, could probe several issues: 



^Filippenko (1988) referred to the events suggested to have hydrogen as Type lb because the Type Ic 
classification (Wheeler & Harkness 1986) was not yet commonly used, but these events are now classified 
Type Ic. 
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e.g., is Ne I A6402 really a candidate to produce the 6200 A absorption, and under what 
circumstances can peaked Ha emission be produced? In the only full non-LTE calculations 
for SNe Ic published so far (Baron et al. 1999), synthetic spectra for models containing no 
hydrogen did not provide good fits to the 6200 A absorption in SN 19941. 

An interesting recent case is that of SN 2005bf (Anupama et al. 2005; Tominaga et al. 
2005; Folatelli et al. 2006), which like SN 1999ex was initially classifed as Type Ic but later 
developed definite He I lines, making it either Type lb or Type Ib/c (or Type Icb?). In each 
of these three papers the 6200 A absorption was attributed to Ha. Strong support for the 
identification was provided by the presence in early spectra of Ca II and Fe II at the same 
high velocity (15,000 km s~^) as the Ha absorption (Folatelli et al. 2006). The very long rise 
time to maximum light of 40 days indicates that the progenitor was massive. In addition, 
SN 2005bf showed evidence for high-velocity polar ejection, reminiscent of SN 1998bw-like 
explosions of massive stars, some of which produce gamma-ray-bursts (GRBs). In view 
of the high ejected mass and the polar ejection, Folatelli et al. suggested that SN 2005bf 
may have been a transition event between ordinary SNe Ib/c and hypernovae. Therefore if 
SN 2005bf ejected hydrogen we should consider the possibility that SN 1998bw-like SNe Ic 
also do. 

In SN 1998bw-like spectra the 5700 A absorption and the 6200 A absorption are assumed 
to be Na I and Si II, respectively. These identifications initially were adopted (Iwamoto 
et al. 1998; Branch 2001) partly because they were assumed to be correct for SN 19941, 
and SN 1998bw-like spectra were and still are regarded to be Doppler-broadened versions 
of SN 19941-like spectra, and partly because Monte Carlo spectrum calculations for bare 
carbon-oxygen cores do provide reasonable fits to SN 1998bw-like spectra (e.g., Mazzali, 
Iwamoto, & Nomoto 2000, Mazzali et al. 2002) with Si II and/or C II producing the 
6200 A absorption. Thus, even if the 6200 A absorption in SN 19941-like SNe Ic is produced 
by Ha, it may still turn out to be produced by Si II in SN 1998bw-likes. Nevertheless, 
the possibility that even in SN 1998bw-likes the correct identification is Ha rather than 
Si II (and perhaps He I rather than Na I) should be considered. Full non-LTE spectrum 
caclulations for comparison with SN 1998bw-like spectra are in their initial stages (E. Baron 
& M. Troxel, in preparation). 



4.2. Implications for Progenitors and Explosion Models 

Ideas concerning the progenitors of SNe lb and SN 19941-like SNe Ic have been reviewed 
by Podsiadlowski (1996), who concluded that most of the progenitors are stars of main- 
sequence masses ^20 solar masses that lose their envelopes by means of binary interactions. 



In standard binary evolutionary scenarios, if all or even just a substantial fraction the helium 
envelope must be lost, all of the hydrogen is lost. For example, Nomoto et al. (1994) favored 
a scenario involving two common-envelope episodes that results in a carbon-oxygen core 
containing little or no helium and no hydrogen. However, as emphasized by Podsiadlowski, 
there are many unresolved questions in binary stellar evolution. If SN 19941-like SNe Ic 
eject hydrogen (and therefore also, of course, helium), some modification of the standard 
evolutionary scenarios, perhaps involving substantial mixing of hydrogen prior to and/or 
during the explosion, will be needed. 

SN 1998bw-like SNe Ic, including those that produce GRBs, are usually regarded to 
come from stars of sufficiently high main-sequence mass that they can lose their envelopes 
by means of single-star winds. Those that produce GRBs may have to be in binaries in order 
to have sufficient angular momentum at the time of core-collapse for the coUapsar model 
(McFadyen & Woosley 1999) to be viable (Fryer & Heger 2005; Petrovic et al. 2005). The 
uncertainities involved in the evolutionary scenarios for SN 1998bw-likes are even greater 
than for SN 19941-likes, and if SN 1998bw-likes eject hydrogen, current scenarios will require 
revision. Particularly interesting in this regard are recent results on chemically homogeneous 
evolution of rapidly rotating massive stars (Yoon & Langer 2005; Woosley & Heger 2006). 
Some of these GRB-candidate models contain hydrogen at the time of core collapse. 

We are grateful to Mario Hamuy, Alex Filippenko, and Tom Matheson for providing 
spectra, to Filippenko, Alicia Soderberg, and the anonymous referee for helpful comments 
on the manuscript, and to Norbert Langer for calling our attention to recent results on quasi- 
homogeneous evolution of rapidly rotating massive stars. This work has been supported by 
NSF grants AST-0204771 and AST-0506028, NASA LTSA grant NNG04GD36G, and NASA 
ADP grant NAG5-3505. 
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Fig. 1. — A day —2 spectrum of the Type Ic SN 19941 is compared with a day —1 spectrum 
of the Type lb or Typelb/c SN 1999ex. The latter is blueshifted by 4000 km s~^ to ahgn the 
absorption features. The flux scale is arbitrary and the spectra have been tilted to facilitate 
the comparison of spectral features. The narrow absorptions near 5900 A in SN 19941 and 
5800 A in SN 1999ex are interstellar Na I in the host galaxies and the narrow absorptions 
near 6700 A and 7400 A in SN 1999ex are telluric. 



-12- 



X 



o 
o 



10 - 



1 1 

19941 day +4 




1 999ex day +4 



3500 



4500 



5500 6500 

WAVELENGTH 



7500 



8500 



Fig. 2.— Like Figure 1 but for day +4 spectra of both SN 19941 and SN 1999ex. The latter 
is blueshifted by 2000 km s~^. 
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Fig. 3. — Like Figure 1 but for a day +10 spectrum of SN 19941 and a day +13 spectrum of 
SN 1999ex. The latter is blueshifted by 2000 km s'K 
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Fig. 4. — The day +4 spectrum of SN 1999ex {solid line) is compared to a Synow synthetic 
spectrum {dashed line) that has Vphot~ ^'^'^'^ ^-^ -^^ 
O I, Mg II, Ca II, and Fe II. 
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Fig. 5. — The day +4 spectrum of SN 19941 {solid line) is compared to a Synow synthetic 
spectrum {dashed line) that has fphot= 11,000 km s^^ and includes hues of hydrogen, He I, 
O I, Mg II, Ca II, Fe II, and Ti II. 
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Fig. 6. — Like Figure 5 except that the synthetic spectrum includes hues of Si II instead of 
hydrogen. 



-17- 



100 



X 

3 



o 



10 - 




4000 



5000 6000 

WAVELENGTH 



7000 



8000 



Fig. 7. — A day +10 spectrum of SN 19941 is compared with spectra of four other SNe Ic 
(from Matheson et ah 2001). The flux scale is arbitrary and the spectra have been tihed to 
facihtate the comparison of spectral features. The spectrum of SN 1990B was obtained on 
day +5. The spectra of SN 1988L, SN 1995F, and SN 1997ei were obtained 11, 10, and 6 
days after discovery, respectively; the dates of B maximum are unknown. 
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Table 1. Synow Fitting Parameters for Figures 4 and 5 
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